Steps in the negative-differential-conductivity regime of a superconductor 
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Current- voltage characteristics were measured in the mixed 
state of YiBa2Cu3C>7-,5 superconducting films in the regime 
where flux flow becomes unstable and the differential conduc- 
tivity dj/dE becomes negative. Under conditions where its 
negative slope is steep, the j(E) curve develops a pronounced 
staircase like pattern. We attribute the steps in j(E) to the 
formation of a dynamical phase consisting of the successive 
nucleation of quantized distortions in the local vortex velocity 
and flux distribution within the moving flux matter. 



In a type II superconductor, a magnetic field H above 
the lower critical value H c \ introduces flux vortices con- 
taining an elementary quantum of flux $ G = h/2e, and 
interactions between the vortices tend to align them into 
a uniform lattice Q]. Extrinsic forces due to impurity 
pinning, thermal fluctuations, or dynamic melting can 
result in a disordered solid or liquid state instead of a 
crystalline lattice 0; however, long-range repulsions be- 
tween the vortices will still enforce a relatively uniform 
density. The system under study consists of a supercon- 
ducting film in a perpendicular applied flux density B a 
along z, with a transport electric current density j and 
electric field E along y in the plane of the film. The trans- 
verse component of E is negligible for this discussion and 
the vortices move with velocity v^, predominantly along 
x. 

A transport current exerts a Lorentz driving force 
Fl = j x $ on the vortices and the motion is op- 
posed by a viscous drag Fd= —V' v 4>i where rj is the co- 
efficient of viscosity. If we assume that pinning forces 
F p are negligible (because Fl 3> F p ) then the steady 
state motion reflects a balance between driving (Fl), 
drag (Fa), and elastic forces (F e ) on each vortex. For 
a perfectly uniform distribution, the net elastic force on 
each vortex vanishes resulting in free flux flow ||] . Then 
j<$> = r)V(j, producing an Ohmic response approximated 
by ij: pf/ Pn « B/H c2 (T). 

A different scenario prevails at ultra high dissipation 
levels and electric fields sufficient to alter the electronic 
distribution function and/or the electronic temperature. 
Here j(E) becomes non-linear and can develop an unsta- 
ble region with negative differential conductivity (NDC) 
(region "C" in Fig. 1). Such unstable behavior and 



NDC has been predicted by Larkin and Ovchinnikov 
(LO) for the regime near T c ||, and has been experi- 
mentally well established ^,0]. In the LO mechanism, 
a non-equilibrium electronic distribution function leads 
to shrinkage of the vortex core by removal of quasipar- 
ticles from its vicinity |^|J^]. In the opposite regime of 
T <C T c we have observed a qualitatively different type 
of instability |||| that seems to result from a tempera- 
ture differential between the electronic system and lattice 
while maintaining an equilibrium-like distribution func- 
tion (because electron-electron scattering is more rapid 
than electron-phonon scattering at low T). In our low- 
temperature instability the vortex expands rather than 
shrinks, and viscous drag is reduced because of a soft- 
ening of gradients of the vortex profile rather than a 
removal of quasiparticles. The standard LO instability 
near T c occurs at values j* and v*—E*jB (see Fig. 1) 
that are B independent, whereas in our low-temperature 
instability j* and v* have a ~ 1/yB dependence. This 
plays an important role in the appearance of a staircase 
at low temperatures but not at the higher temperatures 
of the previously well studied LO phenomenon. Theo- 
retical and experimental details of the low-temperature 
instability are discussed elsewhere p|,|lO|]. 

Fig. 1 suggests a qualitative picture of how the drag 
component (rjv) of the j(E) response might vary with E 
over its entire range. After the first "hump" (regions A, 
B, and C) because of a reduction in rj by the mechanisms 
mentioned above, j(E) will rise again when v$ reaches 
some limiting value such as i>oo ~ £/ta, where £ = 1.5 
nm is the coherence length and ta = h/A 4.7 x lCP 14 s 
is the order-parameter relaxation time. If the vortex den- 
sity is non-uniform, an additional clastic term appears 
in j, namely j = (rjv — F e )/Q . We shall refer to the 
j(E) curve of Fig. 1 as the primitive curve, which applies 
to a hypothetical ensemble of uniformly packed vortices 
moving with identical velocities; an actual j(E) will not 
follow this behavior since, upon encountering a negative 
slope, the flux matter will become unstable and undergo 
a phase transition, such that at any given time there ex- 
ist multiple vortex velocities (resulting in a composite 
response) and a non-uniform vortex density (resulting in 
elastic corrections). In effect some vortices will be travel- 
ling on the second positive slope "E" of Fig. 1 with v<p ~ 
Voc (in regions with reduced vortex density) and the rest 
will be travelling at a reduced velocity v± < v* result- 
ing in restabilization. Below we calculate the expected 
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composite response for one particular flux structure and 
show that it leads to a staircase in a natural way in rough 
agreement with the data. 

The samples are c-axis oriented epitaxial films of 
YiBa 2 Cu 3 7 _5 on (fOO) LaA10 3 substrates with T c 's 
around 90K and of thickness t w 90 nm. Electron-beam 
and optical-projection lithographies, together with wet 
etching in ~ 1 % phosphoric acid, were used to pattern 
bridges of widths w w 4 /jm and lengths I w 90 /an. At 
the end of the lengthy fabrication, each microbridge is 
inspected by a variety of high-resolution optical probes 
for uniformity of width to ensure high reproducibility. 
Altogether ten samples were studied at 12 temperatures 
(1.6, 2.2, 6, 7, 8, 10, 20, 27, 35, 42, 50, 80 K) and at 
11 flux densities (0.1, 0.2, 0.5, 1, 1.5, 2, 10, 11, 13, 13.5, 
14 T). We always observe an instability with steps in 
the NDC region for all temperatures below T c /2 and for 
B values in the 1-14 T range. The electrical transport 
measurements were made with a pulsed constant voltage 
source, preamplifier circuitry, and a digital storage oscil- 
loscope. The pulse rise times are about 100 ns with a 
duty cycle of about 1 ppm, resulting in effective thermal 
resistances of order 1 nK.cm 3 /W. The ability to hold the 
voltage constant across the sample allows investigation of 
the j(E) curve in the NDC region; however, once dj/dE 
becomes negative there will be a jump in the voltage un- 
til the chordal resistance of the sample has risen above 
the source impedance of the voltage source (including 
the resistance of the leads and contacts) allowing a sta- 
ble steady state. Such forbidden gaps in E are visible 
in the last data plot at the onset of NDC. Each j(E) 
curve typically consists of 1000 separate points and each 
point requires averaging over several hundred pulses to 
obtain an adequate signal-to-noise ratio (SNR), so that 
a single curve takes several days to measure. Note that 
the j values in the experiment are an order of magnitude 
lower than the depairing current density jn|. Further 
details about the experimental techniques are discussed 
elsewhere Jl^,|l2|. 

Fig. 2 shows some examples of measured staircase pat- 
terns in the NDC regime. We found such patterns to be 
ubiquitous under most conditions of B a and T (not close 
to T c ). Fig. 2(a) shows the behavior at 11T and 27K. The 
main step features are seen to be reproducible between 
curves measured a week apart in opposite directions of 
changing E. Panel (b) shows 50-K 1-T j(E) curves for 
two different samples (sample X was measured 23 days 
after sample Z). Within the scatter and uncertainly in 
the absolute j values (due to uncertainly in the sample 
widths used to calculate j from I), there apears to be 
some consistency in the main step features. The steep- 
est NDC and most pronounced staircase patterns were 
observed at temperatures around one half of T c and flux 
densities around 1 Tesla (where the intervortex separa- 
tion a ~ y/ & /B becomes comparable to the penetration 
depth). Fig. 2(c) shows 50K curves for sample X at three 
B values over an extended range of E. The solid line fit 



indicates linearity in 1/E over the steep NDC portion. 
Other observed features are a convergence of curves in 
their NDC regions for different i? 's (at fixed T) and a 
growth in horizontal step size with increasing E. 

Let us now consider the nature of the transition in the 
flux matter when the imposed electric field E a exceeds 
E* and enters the NDC region "C" of the primitive curve 
(Fig. 1). Since all the vortices could not be moving uni- 
formly with — E /B (for they would all be unstable) 
the flux matter must reorganize itself by the creation of 
moving defects in the regular vortex structure. The exact 
nature of the defects will depend on the vortex state. In 
a vortex solid there can be pairs of edge dislocations of 
different signs, or defect-interstitial pairs that conserve 
the total flux. If the moving vortex system is a liquid, 
the defect can be a finite magnetic spot or domain with 
a different flux density, since redistribution of a magnetic 
flux in a vortex liquid does not involve significant energy 
barriers (i"3j |. In this Letter we present a simple calcu- 
lation based on the magnetic domain approach, in order 
to provide a simple physical explanation of the effect. 
A detailed and more rigorous approach, which consid- 
ers the distinction among the different possibilities listed 
above, will be presented elsewhere. In the reorganized 
flux structure some vortices at any given time will move 
on the second positive slope "E" with v^, ~ Voo while oth- 
ers move on the first slope "A" with V\ < v* , such that 
the composite macroscopic electric field across a sample 
length averages to E a . The reorganized flux structure 
is influenced by several conditions: (1) The ends of the 
sample are at a fixed potential difference so the average 
electric field along any longitudinal path is close to the 
applied E a ; (2) the current density integrated along any 
longitudinal path is approximately constant due to phase 
coherence (allowing for small variations in the vector po- 
tential); (3) the average flux density in the sample must 
equal the applied value B a , otherwise the demagnetiz- 
ing fields would diverge — this means that if defects are 
formed that have a suppressed local B, the flux density 
B\ in the outside bulk will increase to conserve the total 
flux so as to maintain < B >= Bq; (4) extended low 
B defects will tend to be stationary to reduce circulat- 
ing electric fields around them, which would violate the 
first condition; (5) elastic forces between adjacent vor- 
tices — along with asymmetry in demagnetization (since 
I ^> w) — will tend to discourage relative side-by-side 
motion and promote modulations in v$ and B that tra- 
verse the entire length of the sample; and (6) bulk elastic 
forces within the vortex matter will discourage continu- 
ous growth of defects but instead favour discontinuous 
nucleation of uniformly spaced defects of the smallest 
possible size (~ 1 lattice constant). 

Considering the above guidelines, we estimate below 
the composite j(E) response for one simplified static 
scenario of the defect structure. We would like to em- 
phasize that this is not presented as a comprehensive 
or unique solution to the vortex reorganization problem, 
but only used to provide a simple physical explanation 
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of the observed staircase behavior. Referring to Fig. 1, 
let the applied electric field just exceed the peak value 
E*. The functional form of the primitive j(E) curve in 
the peak region is given by |f||l0| j ~ (a n H c2 / B)E /[\ + 
(E/E*) 2 } with E* = v / 2p n BnA/(H c2 T £ ) and j* = 
\J nAH c2 / (2p n BT e ). Some vortices are forced to move at 
higher velocity in order to restore stability. For the pur- 
pose of obtaining a rough estimate, let's assume that vor- 
tices move at when they cross a longitudinal channel 
with reduced local B and enlarged flux lattice spacing A 
(The actual structure may consist of an array of smaller 
entities that nucleate at the sample edges and migrate 
inward.) The smallest width of this channel, consistent 
with continuity of flux and constancy of longitudinal E, 
is a single lattice spacing A = a(w oc /ui), where a and v\ 
are the lattice spacing and velocity of the "bulk" vortex 
matter outside the channels. As soon as the first defect 
is nucleated, the flux density B\ in the bulk increases 
discontinuously [] from B a to B w/(w — A) and v\ drops 
from E /B to E (w — A) j (B w) . This immediately sta- 
bilizes flux flow in the bulk since E* cx yB increases and 
these vortices drop back on the positive slope "A" (in 
Fig. 1). As the applied E a is increased further nothing 
more happens until v\ reaches the new peak value E* jB\, 
Then a second defect is nucleated and the process keeps 
repeating itself until the defects fill about half the sample 
area, at which point it becomes energetically favourable 
to continuously expand existing defects. The top of each 
step corresponds to j =j*(B\) for each new B\ . Since 
j* oc 1/VB and E* oc \f~B, the staircase shows an over- 
all downward trend with j ~ 1 / E dependence that does 
not depend on B a (causing curves with different applied 
B Q to converge). Both of these behaviors can be seen 
in Fig. 2(c). Fig. 3 shows experimental and numerically 
computed j(E) curves at T=50K and B=l, 1.5, and 2 
Tesla. The two sets of curves show similar qualitative 
trends and have roughly comparable magnitudes. Note: 
The curves computed from the model have no adjustable 
parameters — the values for l/a n ~ p n = 57/xJl-cm and 
H c2 = 120 T come from the literature and E* = 

35 V/cm (at one value of B = 1.5 T) is taken from the 
measurement. 

To summarize, we have investigated the transport re- 
sponse of a superconductor into the regime of negative 
differential conductivity beyond the low-temperature in- 
stability. A qualitatively new behavior was observed in 
the form of steps in the j(E) curve. The observed be- 
havior is consistent with the restabilization of a moving 
vortex distribution by the formation of a dynamical phase 
with distortions in the local flux density and vortex ve- 
locity. While we do not undertake a detailed theoretical 



1 These estimates neglect the relatively small flux 
(BiAZniUoo) through the defects, which, however, is included 
in the numerically calculated curves of Fig. 3(b). 



investigation of the exact nature of the dynamical phase, 
we hope our experimental results will stimulate such fur- 
ther work. Recently other interesting but separate ef- 
fects have been seen in the highly driven vortex state 
by the Huebener group, which include time dependent 
oscillations and hysteretic steps resulting from tunnel- 
ing between neighbouring vortices ffoj ] . The steps in our 
effect are time-independent and non-hysteretic, and are 
analogous to the Gunn effect in semiconductors, where 
electric-charge modulations lead to steps in j(E) in the 
NDC regime. 
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FIG. 1. Schematic diagram for a primitive j{E) curve for 
moving flux vortices distributed with uniform density and 
moving with the same velocity. The net elastic force F e on 
each vortex vanishes and the Lorentz driving force Fl is bal- 
anced by the viscous drag Fd so that j — r]v/<& . Experi- 
mentally only the regime "A" is accessible. Beyond the peak 
at "B" , unstable dynamics destroys the perfect spatial homo- 
geneity of the vortex distribution and leads to non- vanishing 
elastic forces so that j = (r/v — F e )/& . Now vortices travel in 
differing elastic environments and at unequal velocities such 
that each vortex is traveling on either positive slope "A" or 
"E" on the primitive j(E) curve appropriate to its local in- 
stantaneous enviroment. The positive slope "E" arises upon 
reaching some limiting velocity Voo < (,/ta and "F" is entered 
when the sample is driven normal. 
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FIG. 2. (a) Experimental j(E) curves for sample X at 
£? = 11T and T=27 K. The two curve sets were measured 
in increasing (up) and decreasing (down) E to confirm re- 
producibility and absence of hysteresis, (b) Magnified view 
of experimental j(E) curves for two different samples at 50 
K and 1 T (sample X was measured 23 days after sample Z) . 
The step features show some consistency within the scatter of 
the data and uncertainities in the two sample widths used to 
calculate j. (c) 50K low-£? j(E) curves for sample X over an 
extended E range showing overall linearity in 1/E; the solid 
line represents j — 27 + 375/ E.. 
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FIG. 3. (a) Experimental j(E) curves for sample X at 
T=50 K and B= 1, 1.5, and 2 T. In their NDC regimes, 
the curves collapse onto an approximately common behavior, 
where j is roughly linear in 1/E. There are forbidden gaps 
in E just beyond the peak, when the sample's resistance is 
less than the voltage source impedance, as discussed in the 
text, (b) Theoretical j(E) curves for the same T and B's, 
calculated from the simple model based on dynamical flux 
reorganization (the model has no fitting parameters). 
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